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Molecular dynamics simulations are used to study contact between a rigid, nonadhesive, spherical 
tip with radius of order 30nm and a flat elastic substrate covered with a fluid monolayer of adsorbed 
chain molecules. Previous studies of bare surfaces showed that the atomic scale deviations from a 
sphere that are present on any tip constructed from discrete atoms lead to significant deviations 
from continuum theory and dramatic variability in friction forces. Introducing an adsorbed mono- 
layer leads to larger deviations from continuum theory, but decreases the variations between tips 
with different atomic structure. Although the film is fluid, it remains in the contact and behaves 
qualitatively like a thin elastic coating except for certain tips at high loads. Measures of the contact 
area based on the moments or outer limits of the pressure distribution and on counting contacting 
atoms are compared. The number of tip atoms making contact in a time interval grows as a power of 
the interval when the film is present and logarithmically with the interval for bare surfaces. Friction 
is measured by displacing the tip at a constant velocity or pulling the tip with a spring. Both static 
and kinetic friction rise linearly with load at small loads. Transitions in the state of the film lead to 
nonlinear behavior at large loads. The friction is less clearly correlated with contact area than load. 



PACS numbers: 46.55. +d, 62.20.Qp, 81.40. Pq 

I. INTRODUCTION 

Studies of the contact between a spherical tip and a 
flat substrate have played a central role in models of the 
mechanical and frictional response of surfaces. In contin- 
uum mechanics, this geometry can be mapped into con- 
tact between two peaks or asperities on opposing rough 
surfaces as they are pressed together [T]. The response 
of the entire system is then commonly represented as a 
sum over many independent contacts [HEllS], although 
recent studies emphasize the importance of interactions 
between asperities 4; iS] |F| [7]. 

Ideal single asperity contacts are difficult to realize in 
experiments because both the tip and substrate are usu- 
ally rough on scales smaller than the average radius of 
curvature. As a result, most experiments with microme- 
ter and larger tips end up measuring the response of many 
smaller asperities. Notable exceptions are experiments 
on atomically flat mica surfaces in the Surface Force Ap- 
paratus |SJ m Uni- The main limitation of these experi- 
ments is that the contact pressure is typically only tens of 
megapascals, which is orders of magnitude smaller than 
estimates of the pressure between asperities on rough sur- 
faces. 

In the last twenty years there has been great interest 
in the study of single asperities at the nanometer scale 
using the atomic force microscope (AFM) and related 
scanning probes [TTJ [T2 US] • When the AFM is operated 
in a typical contact mode, a tip with characteristic radius 
of 10 to lOOnm is pressed into contact with a surface and 
may be translated to measure the friction. The small 
dimensions of the tip allow more direct control and/or 
measurement of the chemistry and geometry, and reduce 
the range of roughness that can exist on the surfaces. 



AFM experiments provide direct information about 
the normal and lateral forces, but do not reveal the dis- 
tribution of forces within the contact, the area of contact, 
or the mechanisms of deformation and friction. This pre- 
vents them from directly answering long standing ques- 
tions about the role of area and load in determining fric- 
tion. At the macroscopic scale, Amontons's laws of fric- 
tion say that the friction is proportional to load and inde- 
pendent of the area of contact. However a prevailing view 
since the mid 1900's has been that the friction is propor- 
tional to the real area of molecular contact Areai , which is 
much smaller than the nominal or apparent area [U [21 [3] • 
The real area typically grows linearly with load for non- 
adhesive surfaces, and Amontons's laws are recovered if 
there is a constant shear stress Tshear in the contact. For 
adhesive surfaces, A^eai remains finite at zero load, ex- 
plaining why friction is often finite at zero load. 

To try to address the role of Aj-eai at the nanometer 
scale, AFM experiments have been compared to contin- 
uum models. The measured friction can be fit by as- 
suming both that continuum theory describes Areai in 
nanometer scale contacts and that there is a constant 
Tshear [HI El [H] . Howcver, the success of such fits need 
not imply that the underlying assumptions are valid, be- 
cause the functional dependence is quite simple and there 
are many poorly constrained fit parameters. This has 
been demonstrated in subsequent simulations where the 
friction can be fit in the same way, but both assump- 
tions break down [16l [171 HH]- The major cause of the 
discrepancy with continuum theory is the atomic scale 
roughness that is present on any surface composed of 
discrete atoms. This spreads the contact over a larger 
area, much as predicted by continuum calculations for 
rough surfaces pn.[T^. 
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These simulations of single asperity contacts between 
clean surfaces have found a variety of relations between 
Areah load and friction. The friction changes by more 
than two orders of magnitude with the precise atomic 
structure of the tip jl6) , and may scale linearly with area, 
load or neither [HI [171 HI]- Similar strong dependence 
on atomic structure is found for the friction between flat 
surfaces, with friction expected to vanish with increasing 
area unless the surfaces are commensurate, i.e., share a 
common period [20l EJ |22l ESI EH EH ES] . 

Of course experimental surfaces are rarely clean. Any 
surface exposed to air is rapidly coated with adsorbed 
molecules of water, oxygen or small hydrocarbons. Wear 
debris, dust and larger particles may also be present. 
These "third bodies" can have a profound effect on fric- 
tion. Studies of flat surfaces show that they greatly re- 
duce the sensitivity of friction to atomic structure and 
produce a friction force that rises linearly with normal 

load [221 EH EH Eg. 

In this paper we consider the effect of adsorbed mono- 
layers on the mechanical properties and friction of single 
asperity contacts. Even though the adsorbed molecules 
freely diffuse along the surface, they are trapped in a 
glassy state when confined under the tip and remain in 
the contact at very high loads. As for fiat surfaces, the 
adsorbed monolayer decreases the variation of friction 
with atomic structure and yields a force that usually rises 
linearly with normal load. The adsorbed monolayer also 
decreases variations with tip structure in the contact area 
and normal stiffness, but leads to large deviations from 
the results for bare tips. The observed increase in con- 
tact area and decrease in stiffness relative to continuum 
theory are qualitatively consistent with modeling the ad- 
sorbed monolayer as a thin elastic coating. 

The paper is organized as follows. Sec. II describes 
the interactions, geometry and methods used in our sim- 
ulations. The MD results for the pressure distribution, 
contact radius, tip displacement, and friction forces are 
presented in Sec. III. Sec. IV provides a brief summary 
and discussion of the results. 



II. SIMULATION METHODS 

Figure [T] shows a snapshot of the system studied in 
our MD simulations. A spherical tip is pushed, under 
a controlled normal load, into a flat substrate covered 
with an adsorbed fluid layer. The picture only shows a 
small region near the center of contact (about | x ^ in 
the plane of the substrate surface and in the substrate 
depth). The tip radius is i? ~ 30nm, which is comparable 
to radii of tips used in AFM experiments. The geometry 
also approximates contact of peaks on randomly rough 
surfaces. 

In continuum mechanics a contact between two elas- 
tic, frictionless spheres, with radii i?i and i?2, Young's 
moduli El and E2, and Poisson ratios ^1 and V2 respec- 
tively, can be mapped to a contact between a rigid sphere 




FIG. 1: (color online) A tip (top) in contact with a substrate 
(bottom) covered with an adsorbed fluid layer (red). 

of radius R = (i?J~^ -I- ^2^^)^^ and a flat elastic solid of 
modulus E* given by [T] 

1/E* = (1 - uf)/Ei + (1 - iyl)/E2 . (1) 

Previous studies have shown that the mapping remains 
essentially correct at atomic scales [16] . Thus simulations 
can be simplifled to the case of a rigid tip and an elastic 
substrate without loss of generality. 

To approximate a perfectly elastic material, we take a 
face-centered-cubic (fee) crystal as the substrate with a 
(001) surface exposed. Atoms on nearest-neighbor sites 
in the substrate are connected by a harmonic potential 

Vs{r)^hir-d)^, (2) 

where fc is a spring constant, d is the equilibrium distance 
between two nearest neighbors, and r is the inter-atomic 
distance. To be consistent with our previous work, we ex- 
press k and d in terms of the values for particles interact- 
ing with a Lennard- Jones (LJ) potential. The truncated 
and shifted LJ potential is 




for r < Tc, where e and a are the characteristic binding 
energy and diameter, and Vc is the cutoff distance. If 
atoms in the solid interacted with this potential and the 
cutoff only included nearest-neighbors, then one would 
have d — 2^l^a and k approximately ble/a'^. We use 
these values for the spring interaction described in Eq. ^ 
and express all quantities in terms of ct, e and the mass of 
solid atoms m. For example, the unit of time is ^Jma"^ /e. 
The relation between the effective modulus in Eq. ([ij and 
the elastic constants is complicated when the solid is cu- 
bic rather than isotropic fT. The value E* = 63e/cr^ was 
obtained from fits to simulation results for the normal 
displacement in contacts between a bare fee crystal and 
a dense tip where continuum theory is most accurate |17j . 

Periodic boundary conditions are imposed in the x and 
y directions, i.e., in the substrate surface plane. The bot- 
tom layer of substrate atoms is held fixed. This mimics 
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the external support that balances the normal load ap- 
plied to the tip. Continuum solutions assume a semi- 
infinite substrate. To approximate this, all dimensions of 
the substrate should be much bigger than the contact ra- 
dius. The size of the substrate used in this paper is about 
195cr X 195(7 X 190(7, which corresponds to a total of more 
than 6.9 million atoms. The contact radius in our MD 
simulations, as measured by the second moment of the 
pressure distribution, is less than 7% of the substrate di- 
mensions. The quantity most affected by finite substrate 
size is the normal displacement. The magnitude of the 
corrections to the Hertz prediction has been obtained by 
fitting continuum calculations [SHI [511 15^ to simulations 
for bare contacts [17]. The corrections are less than 5% 
and have been included in the Hertz prediction lines for 
normal displacement below. 

We perform a series of simulations using tips with dif- 
ferent atomic scale structures and roughness, but with 
the same radius {R = lOOcr ~ 30nm). The smoothest 
tips are made by bending the (001) surface of an fee 
crystal with nearest-neighbor spacing d' into a spheri- 
cal shape. The friction between bare surfaces depends 
on whether d' /d is rational or irrational. These cases are 
referred to as commensurate and incommensurate, re- 
spectively. The highest friction occurs for commensurate 
surfaces with the same lattice constant and orientation, 
i.e., d' /d =1. In this commensurate case, we find that 
the relative position of the tip on the substrate surface 
plays an important role. Two extreme cases are investi- 
gated here. When tip atoms are directly above interstitial 
positions on the substrate surface, we say the tip is in reg- 
istry since it represents a continuation of the crystalline 
structure of the substrate. On the other hand, when tip 
atoms are directly above atomic equilibrium positions on 
the substrate surface, we say the tip is out of registry. 

Contacting experimental surfaces are rarely in perfect 
registry and alignment. Any mismatch between lattice 
constants (d and d') or misalignment between lattice ori- 
entations leads to incommensurate contact. To model 
such cases, we make an incommensurate tip by bending a 
fee crystal with d' /d = 1.12342. To most closely approx- 
imate a continuum, we also consider a dense tip made by 
bending a (001) face of a fee crystal with d' /d = 0.05 or 
0.1. The two give nearly identical results for contact area 
and normal displacements. Frictional forces are small for 
dense tips and are not reported because such disparate 
lattice constants on opposing surfaces are hard to realize 
with real atoms. 

The tips used in AFM experiments are unlikely to have 
the geometry of bent crystals. To generate more realistic 
tip geometries we take a crystalline or amorphous mate- 
rial and remove all atoms outside a sphere of radius R. 
The crystalline state is chosen to be commensurate with 
d' /d^ 1. The amorphous state is obtained by quenching 
a fluid state with number density l.Ocr"^. While all tips 
considered here deviate from an ideal sphere by less than 
an atomic diameter, the effective roughness is quite dif- 
ferent. The cut crystal has terraced steps on the surface 



and will be referred to as the "stepped" tip in this paper. 
The amorphous surface does not have terraces, but has 
fluctuations in height and density (see Fig. 1 of Ref. [T7]). 

As for the bent commensurate tip, results for the 
stepped tip depend on registry. For most of the results 
presented in this paper the stepped tip is out-of-registry, 
with the atoms in the bottom terrace directly above the 
equilibrium positions of surface substrate atoms. An- 
other case was considered in simulations of friction, be- 
cause sliding along the x axis brings commensurate tips 
between the in and out of registry conflgurations. The 
tip was shifted diagonally by O.bd along both x and y. 
Tip atoms then slide along a line that is centered be- 
tween substrate atoms. This did not change static con- 
tact properties signiflcantly, but lowered the friction by 
a factor of 2. 

Atoms in the tip (t) and substrate (s) interact via 
the LJ potential described in Eq. ([3| with est = 1-Oe, 
CTst = 1.0(7, and Tc — 2^/^(7. The choice of Vc makes 
the interaction purely repulsive. Therefore no adhesion 
is included. Possible effects of adhesion will be explored 
in future studies. Except as noted, the adsorbed film 
prevents direct contact between the tip and substrate. 

The adsorbed monolayer is composed of chain 
molecules containing four atoms. Here "atom" is used 
in a designative sense and each atom may represent a 
monomer of an oligomer chain. Atoms (a) that are not 
connected by a covalent bond interact via the LJ poten- 
tial with Eaa — 0.25e, aaa — 1.0(7, and Vc — 1.8(7. Re- 
ducing the binding energy by a factor of 4 relative to the 
substrate ensures that the molecules melt at a lower tem- 
perature. For adjacent atoms in an adsorbed molecule, 
the LJ potential is truncated at Tc = 2^/^(7, and the co- 
valent bond is modeled by the finite extensible nonlinear 
elastic (FENE) potential 33J 

V.(r) = / -hKRM^ ir/Ro?] if r < i?o, 
^ \ oo if r > i?o- 

Here K — 30eaa/o'^ = 7.5e/(7^ is a spring constant and 
Rq = 1.5(7 is a length scale. Previous studies have shown 
that this coarse-grained model provides a good descrip- 
tion of short-chain hydrocarbon molecules |33j , and it has 
been used successfully to model a wide range of equilib- 
rium and nonequilibrium behavior |34[ 1351 136] . 

The monolayer is confined onto the substrate by an LJ 
potential with esa = 0.4e, asa = 1-2(7, and Vc = I-Sctso = 
2.16(7. The stronger binding energy esa > ^aa gives per- 
fect wetting and, combined with the chain length, leads 
to a very low vapor pressure. All simulations are per- 
formed with the substrate and adsorbate atoms kept at 
a constant temperature T — 0.175e/fcB. This is lower 
than the melting point {T = 0.7e/fcB) of a LJ solid 
whose interactions are comparable to those in the sub- 
strate. However, since eaa is only 0.25e, the tempera- 
ture T — 0.175e/A;B is higher than the glass transition 
temperature (Tg ~ OAeaa/kB ~ O-le/fce) in bulk sys- 
tems [371 ISH] ■ The large asa reduces the periodic poten- 
tial from the substrate that might lock the adsorbed film 
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in an epitaxial state [35]. As a result, the adsorbed film 
is in a fluid state with a high mobility along the substrate 
surface. 

The adsorbates also interact with the tip {t) through 
the LJ potential. The corresponding parameters are 
eta = 0.75e, ata = l.Ocr, and Vc = These pa- 

rameters ensure that the adsorbed film does not wet the 
tip. Thus the formation of a meniscus is avoided. 

The simulations are performed using the Large- 
scale Atomic/Molecular Massively Parallel Simulator 
(LAMMPS) developed at Sandia National Laboratories. 
This classical MD code utilizes spatial decomposition to 
parallelize the computations. Forces are calculated with 
the help of neighbor lists. A velocity- Verlet algorithm 
with a time step dt — O.OOSr is used to integrate the 
equations of motion. Tests with dt — O.OOlr to O.OOTr 
give essentially the same results. To fix the temperature 
T of substrate and adsorbate atoms, a Langevin thermo- 
stat is applied by adding a drag term to their equations 
of motion. The Langevin damping rate F is typically 
0.5r~^. To check that the thermostat did not affect the 
results, we also ran simulations with F = 0.1t~^ and 
damping only on directions perpendicular to the sliding 
velocity in friction measurements. No noticeable change 
was found. 

Although nearest-neighbors in the substrate interact 
with ideal springs, changes in orientation of the springs 
lead to a slight anharmonicity. Instead of expanding 
with increasing T, the equilibrium lattice constant of the 
substrate shrinks by about 0.16% as T rises from to 
0.175e//cB- The periodic boundary conditions prevent 
contraction in the x — y plane, leading to a small negative 
pressure — O.OSe/cr'^) in this plane. The height of the 
substrate shrinks about 0.3% to relieve the tension along 
the z direction. To check that this slight anisotropy did 
not affect our results, we also ran simulations in which 
the solid substrate was allowed to relax to a fully equih- 
brated state with a zero inner pressure before the contact 
occurs. No change was noticed for the results presented 
in this paper. 

A layer of adsorbate molecules is put on top of the 
substrate and the system is allowed to relax into an equi- 
librium state at T = 0.175e/A:B. In Fig. [2] we plot the 
atomic number density p of adsorbates as a function of 
distance h above the mean height of atoms in the top 
layer of the substrate. The film is close to an ideal mono- 
layer with a single major peak at about h ^ l.2a. The 
adsorbate distribution has a long but weak tail at large 
h and a small density peak corresponding to a second 
layer at h ~ 1.9a. Only 2% of the atoms are outside 
the first layer and virtually no evaporation is observed. 
However, the tail in density makes it difficult to define 
the separation between the tip and monolayer covered 
substrate. We will use the minimum in density at 1.65cr 
as a reference height corresponding to the top of the fiuid 
layer. 

After the adsorbed film and the substrate are fully re- 
laxed, a tip is brought into proximity. The normal load 
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FIG. 2: Number density p of adsorbate atoms vs. height 
h above the top layer of substrate atoms before the tip is 
brought into contact. Only points with a nonzero value are 
plotted. Lines are a guide to the eye. The height 1.65a of the 
density minimum (arrow) after the first peak is taken as the 
top of the fluid layer. 

on the tip is controlled and increased by small steps. 
When the contact reaches a steady state at a given load, 
the responses of the adsorbate and substrate are mea- 
sured. The local pressure is determined by finding all 
atoms on the tip or substrate that have a nonzero in- 
teraction with the monolayer or opposing surface. The 
normal component of the force on the atom is then di- 
vided by the area per atom to yield a pressure. We also 
measure the density distribution of the adsorbates, the 
contact area, the tip displacement, and the substrate de- 
formation. Friction is measured by displacing the tip and 
monitoring the lateral force. The tip is moved at a con- 
stant velocity v or attached to a constant velocity stage 
with a spring of stiffness kg. The spring only affects the 
displacement along the pulling direction and no motion 
is allowed along the perpendicular y direction. 



III. RESULTS OF MOLECULAR DYNAMICS 
MEASUREMENTS 

The following subsections compare our simulation re- 
sults to continuum solutions for the geometry studied in 
this paper, a sphere of radius R pushed by a normal load 
N into a ffat elastic substrate. Hertz theory for non- 
adhesive contacts predicts that the region of contact is a 
circle whose radius a is [T] 



The contact area is tto^ and is thus proportional to N^^^. 
The pressure in the contact region takes the form 

, , 2aE* f ^ 
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where r is the distance from the center of the contact 
circle. Hertz theory also predicts that the normal dis- 
placement i5 after initial contact is given by 



S = a^R. 



(7) 



Combining with Eq. [s] this implies that S increases as 

7V2/3. 



A. Distributions of Pressure and Adsorbed 
Molecules 

In Fig. [3] the normal pressure on the substrate sur- 
face is plotted as a function of radial distance r from 
the center of the contact zone under a normal load 
N/{R^E*) = 2.01 X 10-3. The Hertz prediction using 
the bulk modulus of the substrate is plotted as a solid 
line and is the same for all tips. Small dots represent 
the raw data on every substrate atom and circles are the 
average over substrate atoms in an annular region from 
r — 0.5(7 to r -f 0.5cr. Open triangles represent the angle- 
averaged normal pressure on the tip surface. The forces 
on the substrate and tip must balance and from Saint- 
Venant's principal [40^ they are only expected to be re- 
distributed by a distance of order the film thickness. The 
plots show that the two pressure distributions are very 
similar. The only pronounced deviations are near the 
center where fewer atoms contribute to the average. For 
example, only four atoms on the substrate surface are 
located within r < Icr. As a result, the pressure distribu- 
tion is more sensitive to the specific atomic arrangements 
at small radii and shows stronger fluctuations. 
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FIG. 3: Normal pressure distribution in the contact zone for 
different tip geometries: (a) dense, (b) amorphous; (c) in- 
commensurate; (d) commensurate out of registry; (e) com- 
mensurate in registry; (f) stepped. The pressure on each 
substrate atom is shown by a dot in the panels and the av- 
erage over all atoms near a given radius is shown by circles. 
Only the average (triangles) is shown for the pressure on the 
tip atoms. Solid lines indicate the Hertz prediction, which 
is the same in all cases. The normal load is relatively high; 
N/{R^E*) = 2.01 X 10"^ 



The raw pressure data exhibits very strong variations 
with angle. These variations persist even after long av- 
eraging times, but are different for different initial condi- 
tions. Direct observation of the dynamics of the adsorbed 
molecules shows that those in the center of the contact 
region (r/a < 10) are frozen in a glassy state. They un- 
dergo thermal oscillations that are too small to sample 
different pressure distributions in the contact and very 
slow aging that is discussed further below. Similar glass 
transitions have been observed in a wide range of exper- 
iments [HI [HI US] and simulations [HJ on confined 
films. The fiuid outside the contact remains in a mobile 
fiuid state and the dynamics slows down dramatically at 
the outer edge of the contact. 

Previous studies of dry contacts with the same tips 
showed large changes in the distribution of contact pres- 
sure with tip geometry. Fig.[3]shows that tip geometry is 
less important when adsorbed molecules are introduced. 
Tips with very distinct surface configurations lead to very 
similar pressure distributions (Figs. [3^-d). Structural ef- 
fects are only evident for a bent commensurate tip that 
is in registry with the substrate surface (Fig. ^) and a 
tip cut from a crystalline solid (Fig.jsJ'). Moreover, these 
structural effects only become significant at high loads. 
When N/{R^E*) < 2.5 x 10"'', all tips show similar p{r). 
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FIG. 4: Number of monomers per unit area n(r) as a func- 
tion of radial distance r from the central axis (r = 0) of the 
tip for different tip geometries: (a) dense, (b) amorphous; 
(c) incommensurate; (d) commensurate out of registry; (e) 
commensurate in registry; (f) stepped. The normal load 
N/{R'E*) = 2.01 x 10"^ is the same as for Fig.|3]and n(r) is 
normalized by the value before contact no = 0.794a~^. 

To examine how the pressure in the contact affects the 
adsorbed film, the adsorbate surface density n(r) is plot- 
ted in Fig. |4] As in Fig. [3] four of the tips have very sim- 
ilar density distributions (Fig. |4|L-d) even though their 
different atomic structure led to very different behavior 
of bare tips [T5J [T7]. The density in the contact is re- 
duced slightly and shows fluctuations that do not average 
out with time because the film is in a glassy state. The 
density decrease is slightly larger for the amorphous tip, 
perhaps because the greater surface roughness enhances 
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mobility. 

Figure |4|^ shows two dips in n{r) that correspond to 
the edges of terraces on the stepped tip. In continuum 
theory these edges produce a stress singularity, and sub- 
stantial pressure rises are observed for bare tips. The 
adsorbed film smooths these peaks at terrace edges into 
broad features in the pressure plot of Fig. |3f. Note that 
the sequence of terrace sizes is not uniquely determined 
by the tip radius and this would lead to variations in the 
location of these features [T7] . 

Figure ^ shows that the pressure from the commen- 
surate tip in registry expels all adsorbed molecules from 
the central region. The tip atoms interact directly with 
the substrate, leading to a smooth pressure peak in the 
center of the contact (Fig. [3^) that is similar to that for 
a bare surface [ITj . A second pressure peak appears near 
r = 20a where the tip contacts the edge of the adsorbed 
film. For the in-registry tip, atoms from one surface are 
above gaps in the opposing surface. The space available 
to adsorbed molecules is fairly constant and the barriers 
to lateral motion are small. For the out-of-registry case, 
opposing solid atoms are directly above each other and 
adsorbed molecules can be trapped in pronounced chan- 
nels between atoms. We found no expulsion of molecules 
for the out of registry tip even at the highest loads. In- 
registry tips expelled the film for N/{R^E*) > 2.5 x 10"'' 
and there was relatively little dependence on the rate of 
loading because the film was in a glassy state where ther- 
mal activation was negligible compared to mechanical ac- 
tivation. 

All pressure distributions in Fig. |3] exhibit a long tail 
that extends beyond the contact region predicted by 
Hertz theory. Two factors contribute to this. The first 
is that the adsorbed film is compliant and thus can con- 
form to the tip, allowing the pressure at the edge of the 
contact zone to drop to zero more smoothly than in the 
Hertz model. Note that the relatively small density drop 
for most tips indicates that this happens mostly through 
compressing the layer rather than displacing atoms. The 
second factor is thermal fluctuations in the substrate and 
adsorbed film. These are normally neglected in contact 
theories and are discussed in the next subsection. 



B. Measuring Contact Area 

There has been great interest in the relation of single- 
asperity contact area to load and friction [TTl [171 IHl Hll 
SSI iSl SZj. Fig. |3] shows that pressure is distributed 
over a larger area than the Hertz prediction. However, 
the meaning of contact is somewhat ambiguous at the 
atomic scale. Here we will consider three different defini- 
tions of contact radius and area. If we define the contact 
radius as the distance from the contact center at which 
the pressure on the tip surface vanishes (or drops below 
some threshold) , then Fig. [s] indicates that a value big- 
ger than the Hertz prediction will be obtained. We will 
call the contact radius associated with the edge of the 



pressure distribution the outer radius Aq. It is obtained 
from a linear fit to the tail of p{r) on the tip surface. 
Systematic uncertainties are of order a. 

The outer radius is sensitive to the periphery of the 
contact where pressures are small. It is also interesting to 
consider a measure that is most sensitive to high pressure 
regions, such as the second moment of the pressure. If 
the pressure distribution p(r) in the contact zone follows 
the Hertz prediction (Eq. |6]) , then the contact radius a 
satisfies 



2 r^p(r)d^r 
5 Ii^P{r)d'^r 



(8) 



The right hand side of Eq. |8]is readily calculated for the 
actual pressure distributions in our contacts. This gives 
another measure of the contact radius that we will refer 
to as as, where "s" stands for second moment. Green- 
wood and Tripp |19j considered a similar measure in a 
study of the effect of roughness on the contact between 
a sphere and a flat surface using the assumption that as- 
perities could be treated independently. Their radius r* 
corresponded to the ratio of the integral of the pressure 
divided by the integral of p{r)/r. This measure is not 
plotted below because it is less common than the second 
moment and the fluctuations shown in Fig. [3] make r* 
more sensitive to noise than rg. However r* shows the 
same trends as Vg and the deviations are comparable to 
the numerical fluctuations. 

Several groups have attempted to determine the con- 
tact area Ac and the corresponding contact radius Uc 
by counting the number of atoms Nc on one surface 
that are within some cutoff distance of the opposing sur- 
face [TSJ HH 135] . However, this definition is very sensitive 
to the precise atomic structure of the tip and substrate 
and also to the cutoff distance. Often the cutoff is taken 
as the point where the potential becomes repulsive [3H]- 
This criterion for contact correlated most simply with 
friction in a previous study at zero temperature |17j and 
corresponds to the cutoff in the interfacial potential used 
here. 

In our current simulations we found that thermal fluc- 
tuations lead to additional ambiguities in Nc- Fig. [5] 
shows all atoms on an amorphous tip that feel a repul- 
sive force in bare and wet contacts at T = 0.175fcB/e and 
N/{R^E*) = 2.5 X 10"'*. There is clearly a pronounced 
increase in the size of contact when an adsorbed layer 
is present, but there is also a strong dependence on the 
observation time. Different symbols indicate the atoms 
that contacted during time intervals of 50r and 2500r. 
The number of contacting atoms, and thus Ac increases 
by about 25% for the bare contact and by a factor of 2.5 
for the wet contact. 

The variation of Nc with the number of time steps 
t/At is shown in Fig. [6] Values of Nc over a given time 
interval were averaged over many starting times once the 
system had reached steady state. These results are for 
the amorphous tip, but similar results are obtained for 
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FIG. 5: Atoms of an amorphous tip that contact over a time 
interval of 50r (•) or 2500r (□) for a substrate (a) with or (b) 
without an adsorbed film at A'^ = 158e/cr. Large dashed circles 
show the radius ao determined from the pressure distribution. 
This is insensitive to time interval. 

tips with other geometries. For both bare and wet con- 
tcicts the number of atoms in contact is constant at short 
times and then increases monotonically. The length of 
the initial plateau is about a quarter of the period of 
the most rapid vibrations between neighboring atoms, 
i.e., the Einstein modes. This is the minimum time for 
the configuration to change. To confirm this we ana- 
lyzed the autocorrelation function for contact S(t) = 

mp{x,t' + 1)) - {e{p{xm{e{p{x,t')) - {e{p{x)m, 

where p is the local pressure, the Heaviside function 
and angle brackets indicate an average over x and t' . As 
shown for an adsorbed film in the inset of Fig. 16k, S{t) 
crosses zero after about 50 time steps 0.25tP Devi- 
ations from the plateau in grow as S(t) drops and 
the asymptotic long time scaling of Nc sets in at times 
of order 100 to 200 At when S{t) is small. 

The time dependence in Nc comes from several factors. 
For the bare contact the dominant factor is thermal fiuc- 
tuations of the substrate surface. The longer the time 
interval, the more time for a rare event to bring atoms 
at large r into contact. The roughly logarithmic time 
dependence would be consistent with the Boltzmann dis- 
tribution and a linear increase in energy with the height 
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FIG. 6: The number of atoms Nc on an amorphous tip that 
contact the substrate or adsorbed layer over some part of an 
interval containing t/At MD time steps. For both wet (a) 
and bare (b) contacts, the amorphous tip is used. Other tips 
produce similar results. For the wet contact, LogAc oc hogt 
at large t, which means that Nc has a power law dependence 
on t. The slope of the dashed line in (a) is 0.26. For the bare 
contact, Nc increases linearly with logt at large t. The slope 
of the dashed line in (b) is 13.3. The inset in (a) shows the 
autocorrelation function for an atom to remain in contact. 



of surface atoms, but a proper calculation of Nc{t) would 
require a detailed analysis of normal mode energies and 
their effect on contact. For the wet contact, molecular 
rearrangements in the adsorbed film also contribute to 
changes in Nc, leading to more rapid growth with ob- 
servation time: Nc ~ with a ~ 1/4. Slow diffu- 
sion and aging in the center of the contact brings dif- 
ferent tip atoms into contact. At large r, the film is fiuid 
and the cost for thermal fiuctuations that bring adsorbed 
molecules up into contact with the tip is lower than for 
the bare contact. This further broadens the region of 
contacting atoms. 

In the following we consider the contact area associ- 
ated with the initial plateau in Nc, which is just the 
mean number of atoms in contact at any instant in time. 
The corresponding contact area is then Ac = Nc{At)Aat 
where Aat is the area per tip atom. This represents a 
lower bound on the area of contact because of the mono- 
tonic rise in Fig. [6] It would be natural to define the area 
using a longer time interval that averaged out the rapid 
oscillations due to high frequency phonon modes, count- 
ing any atom that contacted over a period of ~ lOOAi. 
Even this short averaging interval increases the area by a 
surprisingly large 50% (25% increase in Cc). This should 
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be kept in mind when comparing to continuum theory, 
but we do not present results for larger time intervals 
because it is hard to identify any unique choice for the 
length of the interval. If the contact region were circular, 
the effective contact radius would be Oc = \/ Ac/tt and 
this quantity is compared to Hertz theory. Note that the 
contacting atoms are in general spread somewhat dilutely 
over a significantly wider region, i.e., ac < Uo- 

Given the clear evolution in the number of contact- 
ing atoms with time interval, one may wonder whether 
there is a corresponding evolution in the contact radii 
measured from the pressure distribution. The answer is 
no. The radius determined from the second moment is 
particularly insensitive to time interval because it is dom- 
inated by the central regions of large pressure. The value 
of flo increases only slightly with time, changing by 0.5a 
or less. The large circles in Fig. [5] show the values of 
determined from the long time interval. Note that some 
atoms outside the circle were in contact during the inter- 
val. However, the magnitude and duration of the forces 
during these contacts were so small that they did not 
contribute significantly to the average pressure and thus 
Oq. Conversely, while many of the atoms inside Oq do not 
contact during the short time interval, we find that there 
are always a few atoms with r w Oq that make strong 
enough contacts to produce a significant pressure. The 
identity of the atoms that make this contribution varies 
from one time interval to the next, but the measured Uq 
is nearly unchanged. 



C. Variation of Contact Area with Load 

Our results for the load dependence of different defini- 
tions of the contact radius are shown in Fig. [7] Values for 
the bare amorphous tip and the Hertz prediction (solid 
line) are shown for comparison. In most cases the radii 
rise linearly with N^^^ and the slope is close to the Hertz 
prediction. However, there are large shifts between dif- 
ferent measures of contact radius. 

The radius Qq where the pressure goes to zero (Fig. [7|l) 
shows the largest deviations from Hertz theory. A sig- 
nificant increase over Hertz theory is observed for bare 
tips [16l [49] . Adding an adsorbed film spreads the pres- 
sure over a much larger area, as already seen in Fig. [5] 
Indeed the value of Uo at the smallest load is bigger than 
the value predicted by Hertz theory for the largest load 
we studied. Note that at this radius, a — 0.14i?, the tip 
surface is only O.Oli? above the lowest point on the tip. 
For R — lOOcr this height is comparable to the width 
of a molecular layer and it is not surprising that the 
film can spread contact over this distance. This effect is 
likely to be significant for most AFM tips, but the broad- 
ening would be insignificant for the much larger radii 
(i? ~ 10mm) in Surface Force Apparatus measurements. 
The shift in Qq for bare tips can also be attributed to a 
deviation from hard sphere behavior. However, instead of 
having an adsorbed layer that spreads interactions over a 




(3N/4E*R2)i« 

FIG. 7: Variation of contact radius a with normal load A'" for 
different tip geometries: dense (x), amorphous (0)i incom- 
mensurate (□), commensurate out of registry (A), commen- 
surate in registry (0), and stepped (v)- Results from Hertz 
theory (solid lines) and for a bare amorphous tip (-I-) at the 
same temperature {T = 0.175e/fcs) are shown for compari- 
son. Three methods are used to measure the contact radius: 
(a) Ho from the edge of p{r); (b) as from the second moment 
of p(r); (c) ac from the mean number of contacting atoms at 
each MD time step. 

height of order a, the broadening comes from the shorter 
lengths associated with the range of repulsive interactions 
{r-^ 0.1a) and amplitude of surface roughness [T6tl49]. 

For the stepped tip, rises in discrete steps as N 
increases. As shown in Fig. |7|l, it is almost constant for 
the lower three loads, jumps to another plateau for the 
next two loads, and has another discontinuous increase at 
the largest load. This is due to the fact that the surface 
of this tip has terraced steps. Each time that the next 
terrace comes into contact with the adsorbed film, Oq 
increases to the outer radius of that terrace. Note that 
the sequence of terrace sizes is not uniquely determined 
by the tip radius and thus the load where jumps occur 
would vary with realization. 

The pressure in the outer regions of the contact is small 
(Fig. [4]) and does not contribute significantly to the load 
on the tip, particularly at large N. Fig. ^ shows that 
the contact radius determined from the second moment 
of the pressure is much closer to the Hertz prediction. 
The results at the largest loads are generally parallel to 
the Hertz prediction and shifted upwards by only one 
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or two molecular diameters. The only exception is the 
commensurate tip in registry. As noted above, it expels 
the film at the highest three loads and there is a sec- 
ond peak in the pressure at the edge of the expelled film 
(Fig. |4^) that increases as- The difference between bare 
and wet contacts for amorphous tips is small, indicating 
that most of the contacting atoms in the outer region 
of Fig. [SJd carry relatively little load. Similar results are 
obtained when comparing bare and wet contacts of other 
tips. 

At the lowest loads the value of saturates or actually 
increases with decreasing load. This counterintuitive be- 
havior results from the very low pressures in the contact. 
The lowest load is 5.87e/tT, corresponding to a pressure 
of only 0.04e/cr^ when spread over the observed radius of 
order 7a. The thermal contribution to the pressure in a 
system of density n is nk-gT. Thus the lowest load could 
be balanced by an ideal gas with density of only 0.22cr~'^ 
for the temperature in our simulations. Hertz theory ig- 
nores thermal fluctuations and can not be expected to 
apply in this limit. The pressures corresponding to the 
second and third loads are about 8 and 27 times higher, 
respectively. The thermal contribution to the stress be- 
comes increasingly irrelevant as the load increases and 
the Hertz prediction becomes more accurate. Note that 
at the lowest load we also find that the adsorbed layer is 
not in a glassy state. 

It is interesting to consider how the loads considered 
here compare to typical experiments. If we take E* = 
lOOGPa and R = 30nm, the loads considered here would 
range from about 1 to 200nN. This is comparable to 
the range of experimental loads. However, experimental 
systems are often adhesive. Even weak van der Waals 
interactions lead to forces per unit area that are of order 
k^T /a^ or larger and this may limit the role of thermal 
fiuctuations. Studies of this are underway. 

Figure [TJ: shows results for the radius Oc, correspond- 
ing to the area of atoms in contact at any instant. As 
for as, the increase in Oc is roughly parallel to the Hertz 
prediction, but ac is generally shifted to lower radii. As 
noted above, ac is expected to be smaller than other mea- 
sures. Including all atoms that contacted on the time 
scale of high frequency vibrations would increase Oc by 
about 25% and make it comparable to at large loads 
for most tips. 

The results for the dense tip indicate another difficulty 
that arises in defining A,.- Even for ideal flat surfaces, 
only a small fraction (~ 1/6) of the atoms on the dense 
tip are close enough to contact the opposing surface. 
Atoms that lie between atoms on the opposing surface 
are too far away to contact. Normalizing by the fraction 
of atoms that are in contact for two flat surfaces gives 
values of a^ that are quite close to those for other tips. 
Similar effects are likely to arise on surfaces with a range 
of surface species, where some species may be screened by 
others with greater height or larger interaction lengths. 
The relative contribution of different species to contact 
area is also unclear in this case. 



Note that the results on Oc for the in-registry commen- 
surate tip rise slightly above the Hertz prediction at the 
three largest loads. This is because the tip breaks the ad- 
sorbed fllm and makes direct contact with the substrate. 
If we only count the number of tip atoms in that direct 
contact, then ac would be smaller than the Hertz predic- 
tion. However, including contacts with the film at the 
edge of the contact zone gives a larger ac- 

D. Normal Stiffness of the Contact 

Experiments can not measure the pressure distribution 
or contact radius directly, but the normal displacement 
can be determined if the device is sufficiently stiff [12]. 
Hertz theory assumes hard sphere interactions and the 
tip displacement is defined as the indentation 5 relative 
to the first contact. In experiments and simulations the 
interactions always have a finite range that leads to am- 
biguity in the zero of 5. In our simulations, the long tail 
in the density shown in Fig. [2] leads to weak repulsions at 
large separations. As noted in the discussion of this fig- 
ure, we took the density minimum at 1.65(j as the outer 
edge of the first monolayer. The tip would then contact 
this layer when height of the lowest tip atom was higher 
than this outer edge by the cutoff radius 2^/^cr. We use 
this height as our zero for 5. 
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FIG. 8: The normal displacement of the tip & (filled symbols) 
and of the substrate &s (open symbols) vs. load for the in- 
dicated tip geometries. The solid line is the Hertz prediction 
with E* — 63e/cr^, R — lOOcr and a correction for finite sub- 
strate size [17]. The dotted line shows the uncorrected Hertz 
prediction and the dashed Hue is the corrected result with a 
constant offset. 

The filled symbols in Fig. |8]show the normal displace- 
ment of the tip b for different tip geometries. The mean 
tip height was determined by averaging over at least 
2000t at each fixed load. Uncertainties due to thermal 
fiuctuations were less than O.OScr, and greatest at low 
loads. The results are plotted against A^^/-^ so that the 
Hertz prediction for an infinite substrate is a straight 
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line. This is indicated by a dotted line. A solid line 
shows the continuum prediction with a correction for fi- 
nite substrate depth. The correction is less than 5% and 
smaller than the symbol size. 

For all tips, the value of S rises much more rapidly 
than predicted at low loads. At higher loads, the rate of 
increase slows and results for most tips can be fit by the 
Hertz prediction with a constant offset. In experiments 
the modulus of a substrate is often obtained from the 
slope of plots of S against N^^^. Based on the results 
shown in Fig. [8] such fits should give reasonably accurate 
moduli 10%) if extended to sufficiently high load. 

Similar load-displacement curves are obtained in cal- 
culations for substrates that are covered by a thin layer 
that is more compliant [311 132] • When the contact radius 
is much smaller than the layer thickness, indentation is 
described by the lower modulus of the layer and 5 rises 
more rapidly. When the contact radius is much larger 
than the thickness, the effective compliance is that of 
the stiffer substrate. Plots like Fig. |8] should show two 
straight lines at low and high loads with slope changing 
by the ratio of E*'^^^. The change of slope in our calcu- 
lations by almost a factor of three would imply an order 
of magnitude increase in modulus. Indeed, as discussed 
above, the film may be more fluid than solid at the lowest 
load. 

To isolate the amount of displacement accommodated 
by the film, we also evaluated the mean normal dis- 
placement of atoms in the top layer of the substrate 5s ■ 
Here the origin can unambiguously be taken as the mean 
height in the absence of any load. Moreover, Eq. [7] cor- 
responds exactly to the expected substrate deformation 
for a rigid tip. 

As shown in Fig. [8] 5s is consistently below the Hertz 
prediction for all tips. The reason is that the adsorbed 
film distributes the force from the tip over a wider area. 
This increases the effective stiffness of the substrate and 
lowers the displacement. The broadening of the pressure 
becomes less important as the contact grows and 5s is 
parallel to the Hertz prediction at large loads. The in- 
registry commensurate tip has a very large 5 because 
the adsorbed film is pushed out. However, it produces 
a slightly smaller 5s than those of other tips. This is 
because part of the load is spread to the fluid around the 
tip and spreading the load over a larger area makes the 
system stiffer. 



E. Time Dependence of Friction 

Measurements of friction are known to be affected by 
the mechanical properties of the measurement device, 
whether at macroscopic [5] or nanometer [501 Ell El] 
scales. To illustrate the range of possible friction forces, 
we contrast the friction measured by stiff and compliant 
systems. In the stiff case, the tip is displaced at a con- 
stant velocity v along the x direction and the only compli- 
ance is provided by the substrate and the interface with 
the tip. In the compliant case, the tip is connected to one 
end of a spring of stiffness ks and the other end is dis- 
placed at constant velocity. The value of ks for each tip 
is chosen so that the spring is more compliant than either 
the substrate or contact. To limit the effects of inertia 
that are described below, a damping force is added so 
that the spring is critically damped. The damping does 
not contribute directly to the reported frictional forces. 
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FIG. 9: Variation of friction force with time for (a) out-of- 
registry commensurate and (b) amorphous tips driven with 
constant tip velocity (solid lines) and through a compliant 
spring (dashed lines). Since the velocity is O.Olcr/r, a time 
interval of 150r corresponds to a displacement by the period 
ci = 1.5(T of the substrate along the direction of motion. The 
normal load is A'' = 1267.2e/(j and the spring compliance is 
ks = 40e/cr^ 



The total change in 5 is larger than that in 5s- The dif- 
ference must be accommodated by compressing the ad- 
sorbed film and the interfaces on either side. For the 
cases studied this is typically less than 0.5cr. This to- 
tal change should be insensitive to the radius of the tip. 
Thus normalized plots like Fig. [8| will show larger changes 
for smaller R and smaller changes for larger tip radii. 



Figure |9]shows the the time dependent friction force on 
out-of-registry commensurate and amorphous tips driven 
by stiff or compliant systems at w = O.OIct/t. The com- 
mensurate tip exhibits simple, nearly periodic motion 
that is often seen in experiments jllj . Due to the periodic 
structure of the tip, the system can find an equivalent en- 
ergy minimum after a displacement by a lattice constant 
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d. The amorphous tip is disordered and its motion is less 
regular. 

The maximum friction force corresponds to the static 
friction Fg that is needed to initiate steady sliding. For 
both tips the heights of force maxima are relatively in- 
sensitive to the compliance of the system. However the 
heights do vary with position, particularly for the amor- 
phous tip. These fluctuations reflect changes in the con- 
formation of the adsorbed layer and tend to increase in 
magnitude with velocity [53]. The kinetic friction F^. 
corresponds to the time average of the friction force and 
varies strongly with compliance. Stiffer systems have a 
lower Fk because the friction minima decrease with in- 
creasing stiffness and may even become negative. 

To understand these trends it is useful to consider the 
simple Prandtl-Tomlinson model |50|, I54[ 155] for motion 
of a single degree of freedom, the rigid tip, in a potential 
with period d. The equation of motion is 

2nx 

mx — ks(xo — x) — Fq sin — Tx , (9) 

where m is the mass, x is the location of the tip, dots 
indicate time derivatives, kg is the stiffness of a pulling 
spring whose other end is at position xq , and F is a damp- 
ing coefficient. The static friction is just the maximum 
force exerted by the periodic potential, Fq. The kinetic 
friction exhibits two different types of behavior depend- 
ing on whether kg is smaller or larger than a critical value 
kc = 2TTFo/d. 

For kg > kc there is a single static solution to Eq. |9]for 
each xq. As Xq = vt increases, the tip moves smoothly 
with X V and the friction from the periodic potential 
varies from positive to negative values. The resulting 
time average force is small and Fk goes to zero with de- 
creasing V. This type of behavior is exhibited by the com- 
mensurate tip when the only compliance comes from the 
substrate. Note that the negative minima in the force 
are nearly equal and opposite to the positive maxima. 
The time averaged kinetic friction Fk = 34e/a is less 
than 10% of the static friction Fg — 368e/a. To conflrm 
that kg > kc for this system, we note that the values of 
Fg and d imply kc — 1462e/CT^. The compliance of the 
substrate is comparable to the prediction of continuum 
theory kg = 8Ga, where G = 18.3e/a^ is the shear mod- 
ulus and a the contact radius. Since all measures of a 
discussed above are larger than the continuum predic- 
tion of 11.5(7 for this load, kg > 1684e/(T^ > kc. The 
same condition holds for all other loads because Fg de- 
creases more rapidly with load than the contact radius 
does. 

When kg < kc there are multiple static solutions to 
Eq. |9] for a given xq. The tip remains stuck in one un- 
til it becomes unstable and then jumps rapidly to an- 
other static solution. The force is near Fg before the 
jump, and only drops by of order kgd since the distance 
between stable states is of order the period. As a re- 
sult, the time average kinetic friction approaches Fg as 
fcs — > 0. For the commensurate tip with a compliant 



spring, kg = 40e/cr2 « kc and Fk is almost 90% of Fg. 
Studies with weaker and stronger springs confirmed that 
Fk I Fg varies from nearly zero at kg > kc to almost unity 
in the limit kg « kc. 

The time dependence of the compliant system illus- 
trates another effect that can change the measured ki- 
netic friction. In Fig. [9j^a) most of the drops in force 
have the same size, but two later drops are larger by a 
factor of three (t ^ 1400t) and two {t ~ 2200t), respec- 
tively. In these cases the inertia of the tip allowed it to 
jump forward by more than one period. These multiple 
slip events are also observed in some AFM experiments 
and have been modeled using Eq.[9][FIllIl]. They lead 
to a substantial reduction in Fk by reducing the minima 
in the force. 

The frequency and magnitude of multiple slip events 
vary with load as well as compliance and tip iner- 
tia [511 US- While this may be important in certain 
experimental systems, it makes it difficult to report a 
single number for the load dependent kinetic friction on 
our model tips. To minimize the associated ambiguity, 
we will focus on cases where multiple slips have been 
eliminated by critical damping of the spring and use of a 
sufficiently stiff pulling spring for each tip. 

The time dependence of the friction force on the 
stepped and incommensurate tips is qualitatively simi- 
lar to that of the commensurate tip. In both cases, the 
tip can flnd an equivalent minimum if it moves forward 
by a lattice constant of the tip while the adsorbed fllm 
remains unchanged. Note that translating the incom- 
mensurate tip produces a different registry with the sub- 
strate, but this does not matter if the adsorbed film has 
already adapted to the periodicity of the tip. 

In contrast, the time dependent friction on the amor- 
phous tip is not simply described by Eq. [9] Because the 
structure of the tip is not periodic, the tip does not find 
an equivalent potential energy minimum after moving by 
any distance. Periodic motion would be possible if the 
adsorbed molecules slid rigidly over the substrate with 
the tip, but this is not observed. At low loads, sliding is 
localized at the interface between the tip and film. The 
film is damaged at higher loads, as discussed in the next 
subsection. While the tip drags some adsorbed molecules 
over the substrate, they do not translate rigidly with the 
tip and thus do not lock into an equivalent energy mini- 
mum after moving by a period. 

The lack of periodicity has a direct impact on the fric- 
tion traces for amorphous tips. Even for the stiff system, 
where kg > kc = 420e/cr^, the friction force remains pos- 
itive as the tip jumps forward. As a result, Fk is close to 
Fg for any degree of compliance. The minima in friction 
are typically about half the maxima rather than being 
equal and opposite in magnitude as predicted by Eq. |9] 
The negative forces for the stiff commensurate system 
come from the intervals where the tip is being pulled for- 
ward as it drops into a free energy minimum. There is no 
preformed minimum for the amorphous tip to drop into 
and thus no negative force. 
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For the stiff amorphous system the force shows oscil- 
lations that are not periodic but have a typical spacing 
that is of order the atomic diameter. This length scale 
enters simply because it is the typical distance that the 
tip must move to interact with different molecules af- 
ter leaving a low free energy state. For the compliant 
system the motion is even less regular. There are some 
intervals where the force is fairly constant, and the tip 
moves over the substrate at a nearly constant velocity. 
At other times the tip gets stuck in a local free energy 
minimum. The force then rises linearly with time un- 
til the minimum becomes unstable and the tip jumps 
forward. Similar erratic stick-slip has been observed in 
previous simulations and was found to be sensitive to the 
inertia of the tip [53J . The simulations reported here are 
overdamped, which decreases the role of inertia. 



F. Load Dependence of Friction 

The static and kinetic friction were obtained from fric- 
tion traces like those shown in Fig. |9] The kinetic fric- 
tion is just the time average of the force after the system 
reaches a steady state. The static friction was obtained 
by identifying the peak forces during each interval where 
the surfaces locked together. We found the peaks were 
easiest to identify for the compliant systems and present 
those results below. Values obtained from stiff systems 
were the same within the statistical uncertainty, which 
was less than 10%. Experiments, previous simulations, 
and simple models suggest that the static friction should 
grow with the time the surfaces are stuck together be- 
cause of aging in the adsorbed film |3S1 |3SJ [571 EH] . This 
dependence is typically logarithmic in contact time and 
difficult to separate from statistical fluctuations. It is not 
considered here. 

The static and kinetic friction were evaluated over a 
sliding distance of at least 20(t. At low loads the fric- 
tion varied little over this interval. At the highest loads 
some tips induced transitions in the state of the film with 
sliding distance that changed the friction force. 

The evolution of film structure at = 1267.2e/(T is 
At this load, there is a tendency 



illustrated in Fig. 10 



for a hole to form behind the tip and for molecules to pile 
up in front of the tip. The size of these effects increases 
with the roughness of the tip. The smallest changes are 
for the bent incommensurate and out-of-registry com- 
mensurate tips. Both slide on top of the film at all loads, 
with almost all of the sliding velocity accommodated at 
the tip/film interface. In contrast, the stepped tip com- 
pletely expels film molecules from under the first ter- 
race for N > 375e/(T. Molecules remain trapped under 
the second terrace and tend to move with the tip. For 
N — 1267.2e/CT, even the adsorbed molecules under the 
second terrace start to be squeezed out. This allows the 
first terrace to contact the substrate directly for the first 
time and leads to a rapid rise in friction. 

The amorphous tip slides on top of the film at low 
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FIG. 10: (color online) Snapshots showing the location of 
all film atoms (dots) and those tip atoms that contact the 
film (open red circles) or substrate (filled blue circles) over 
an interval of 0.5r before (left) and after (right) the tip is 
pulled over the substrate. The tip geometries are: out-of- 
registry commensurate (a and b); incommensurate (c and d); 
amorphous (e and f); and stepped (g and h). The vertical 
height of the panels is 50a and the horizontal dimension is 
the same for the left panels and 85a for the right panels. 
Large circles show the contact radii determined in Fig. [7] with 
the largest corresponding to Uo and the smallest to a^. The 
only place where tip atoms directly contact the substrate is 
within the inner circle in panel (g). 



loads, but the mode of sliding changes for N > 733e/cr. 
Some of the molecules become trapped in the roughness 
of the tip and are dragged along with the tip. The slid- 
ing velocity is accommodated between these molecules 
and the substrate. Because the molecules are trapped in 
positions related to the roughness of the disordered tip, 
they do not lock well with the substrate and the friction 
rises less rapidly with N than at lower loads. The tip 
and entrapped film molecules push the rest of the film 
aside. Molecules pile up in front and a bare patch is left 
behind. This bare patch fills in fairly rapidly and the re- 
gion where the tip started has nearly healed by the end 
of the run. Such healing of the lubricant layer could be 
very important to function as a boundary lubricant. 

A different type of time dependence was observed for 
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the in-registry commensurate tip at high loads {N > 
375e/(T) and velocities. For v > O.Olcr/r the tip starts 
in direct contact with the substrate and the first friction 
peak is about twice as large as for the out-of-registry tip. 
Once sliding starts, the tip lifts up onto the film. The fric- 
tion drops sharply and then rises gradually, with the final 
friction being very close to that for the out-of-registry 
commensurate tip. For this reason we do not show sepa- 
rate friction results for the in-registry commensurate tip. 
When the velocity was decreased to u = 0.002cr/T the tip 
remained in contact with the substrate during the slid- 
ing. In this case the tip plows through the film and the 
force peaks remain roughly twice the height of those at 
V = O.Olcr/r. 

For both the in and out of registry bent commensurate 
tips, the sliding path takes atoms on the tip directly over 
atoms in the substrate. This is why the frictional forces 
are identical as long as the film is either present or absent 
for both tips. Previous work has shown that the friction 
on commensurate tips changes as the line of motion is 
displaced perpendicularly to the sliding direction [59l [Ml 
I5T| . To illustrate the effect of perpendicular shifts, the 
stepped tip was displaced diagonally by 0.5d along x and 
y so the tip atoms move along a line that is centered 
between substrate atoms. Before displacing the friction 
on the stepped tip is close to that on bent tips. After 
displacing the friction is reduced by a factor of about 
two, with the precise ratio depending on load. Similar 
reductions are observed when the bent commensurate tip 
is offset in the same way. The stepped tip was also made 
incommensurate by rotating it around the z-axis by 20° . 
This produced an even larger drop in friction. 

Figure [TT] shows the load dependence of the static fric- 
tion (Fg) on different tips and the kinetic friction (Fj.) 
obtained for compliant and stiff systems. While the fric- 
tion varies with tip geometry, the deviations are much 
smaller than those observed for bare tips. For the same 
tips, the friction in bare contacts varied by two orders of 
magnitude and showed different functional dependencies 
on load. Inserting the adsorbed layer has greatly reduced 
this sensitivity to tip geometry, and the variation is even 
smaller if the special case of commensurate tips moved 
directly over substrate atoms is removed. In addition, 
the tips all exhibit a roughly linear rise in friction with 
load, although there are transitions in behavior for some 
tip geometries. For example the friction on the stepped 
tip increases more rapidly when the film is expelled from 
the central terrace {N > 375e/f7), while the friction on 
the amorphous tip increases less rapidly with load when 
the film is damaged [N > 733e/(T). 

In these and previous simulations the largest static fric- 
tion is observed for the commensurate case where the 
opposing surfaces have the same period and orientation. 
This case also typically yields a linear relation between Fg 
and load. The simplest explanation is that the surfaces 
are effectively in hard-sphere contact and the tip must be 
lifted up over atoms in the layer below before it can slide 
forward [55^. The friction provides the force needed to 
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FIG. 11: Static {Fs) and kinetic (Ffc) friction vs. load (iV) 
for the indicated tip geometries: amorphous (0)i incommen- 
surate (□), out-of-registry commensurate (A), and stepped 
(v)- The static friction in (a) and the kinetic friction in (b) 
are measured by dragging the tip along the x-direction with 
a spring at a constant velocity O.Olcr/r. The kinetic friction 
in (c) is measured by displacing the tip uniformly along the 
a;-direction at a velocity 0.002(T/r. Straight lines are linear 
fits to the friction. Except in (c) , they go through the origin. 



move the tip up the ramp formed by these atoms and the 
slope of the ramp gives dFs/dN . The slope is largest for 
the commensurate case because the effect of all atoms 
adds coherently. The slope is also larger when the tip 
atoms pass directly over the substrate atoms rather than 
passing between them as for the stepped tip simulations. 
For incommensurate and amorphous surfaces, the forces 
average to zero for bare, infinite surfaces. This cancella- 
tion is prevented by the adsorbed film because it adopts 
a glassy state that adapts to the structure of both tip 
and substrate [57]. However, the commensurate case still 
proves most easily locked together. 

The kinetic friction for compliant systems is very sim- 
ilar to Fg. The stiff amorphous results are also com- 
parable to Fs for the reasons discussed in the previous 
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section. The friction force on commensurate systems is 
substantially reduced in the stiff system. The friction 
is nearly zero at small loads and then rises linearly at 
large loads. This behavior can be understood from the 
Prandtl-Tomlinson model. At low loads, kg > kc and the 
kinetic friction vanishes with decreasing velocity. As the 
load increases, the potential and k^ increase. At suffi- 
ciently large load, kg may be smaller than kc leading to 
a kinetic friction that follows the linear rise in Fg. As 
expected from Fig. [9] the stiff commensurate system has 
kg > kc for all loads. Even though the commensurate tip 
has the largest static friction of any tip, Fk is the smallest 
and goes to zero with decreasing velocity. 

There has been great interest in understanding the con- 
nection between the friction in molecular scale asperities 
and macroscopic measurements of friction. Amontons's 
laws state that macroscopic friction is proportional to 
load and independent of the nominal area of the contact- 
ing surfaces. However, macroscopic surfaces are generally 
rough and, as noted above, the real area of molecular 
contact may be much smaller than the nominal area. In 
many cases, Areai is expected to be proportional to load 
but is difficult to measure, making it unclear whether 
the area or load is controlling the friction force. There 
are also many exceptions to Amontons's laws. For exam- 
ple, adhesion between surfaces leads to friction at zero 
or negative loads and grows in importance as dimensions 
shrink. Repulsive tip interactions were used in our sim- 
ulations to eliminate adhesive effects. 



thors reported the static friction while Mo et al. reported 
kinetic friction and went to much smaller dimensionless 
loads. 



The variation of friction with load in Fig. 11 is roughly 
linear and thus consistent with Amontons's laws. As 
noted above, deviations from linearity tend to coincide 
with changes in the structure of film and location of slid- 
ing. Linear fits to the kinetic friction reach zero fric- 
tion at a positive load, but this shift depends on the 
pulling spring and can be understood from the Prandtl- 
Tomlinson model (Eq. [9]) . 

Studies of bare tips by Luan and Robbins showed that 
the friction did not always increase linearly with load over 
the same range of loads studied here [T7]. In particular, 
the static friction on the amorphous and incommensurate 
tips increased sublinearly and the stepped tip showed dis- 
continuous jumps in friction. When a film is present, the 
friction on the amorphous tip only rises sublinearly at 
the highest load and this is directly attributable to a 
change in sliding mechanism that one would not expect 
to describe simply in terms of the change in load or area. 

Mo et al. have recently considered the friction between 
diamond and a hydrogen-terminated amorphous carbon 
tip, using realistic interaction potentials [18^. They con- 
cluded that both the friction and the area obtained from 
counting atoms, Ac, were proportional to load and con- 
trasted this with the results of Luan and Robbins. The 
plots shown by Mo et al. actually show that neither fric- 
tion nor area is strictly proportional to load. Rather both 
are zero up to a finite load and then rise linearly. While 
this linear relation is still inconsistent with the sublinear 
relation obtained by Luan and Robbins, the latter au- 
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FIG. 12: Friction on a bare amorphous tip as a function of 
dimensionless load at v = O.Ola/r and ksT/e = 10~*. The 
static friction (open symbols) and kinetic friction (filled sym- 
bols) are shown for constant tip velocity (squares) and springs 
with stiffness 40e/a^ (downward triangles), 2Qe/a^ (upward 
triangles), and lOe/cr^ (circles). 

Figure [T2| contrasts the static and kinetic friction on a 
bare amorphous tip over the range of dimensionless loads 
considered by Mo et al.. The static friction is strongly 
sublinear as noted by Luan and Robbins. The kinetic 
friction is consistent with the results of Mo et al. . As ex- 
pected from the Prandtl-Tomlinson model, the friction is 
nearly zero at low loads where the driving system is stiff 
compared to the substrate potential. At larger loads the 
friction rises roughly linearly with load. The crossover 
between these behaviors moves to larger loads as the 
stiffness of the pulling spring increases. Note that the 
sensitivity to stiffness is lower when an adsorbed film is 
present because the potential energy is not a simple func- 
tion of tip position as assumed in the Prandtl-Tomlinson 
model. 

The results in Fig. [12] explain the apparent discrep- 
ancy between previous results. Luan and Robbins did 
not report Fk because of the strong dependence on stiff- 
ness. Mo et al. did not report Fg or vary the system 
compliance by adding a spring. However, the tip, sub- 
strate and hydrogen bonds terminating the tip all act as 
compliant elements. The importance of grafted hydrogen 
and hydrocarbons in locking surfaces together and dissi- 
pating energy has been shown in a number of studies by 
Harrison and coworkers [511 . 

Both previous studies of bare tips found that the fric- 
tion rose roughly linearly with area as measured either 
by tto [ISJ EZ] or Oc [H]- However, the friction did not 
necessarily got to zero at zero area, and in many cases 
the outer region contributed little to the friction. Fig. [13] 
shows a plot of friction vs. for tips sliding on a mono- 
layer. The behavior at large forces could be fit to a line 
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squared for the indicated tip geometries: amorphous (0)i 
incommensurate (□), out-of-registry commensurate (A), in- 
registry commensurate (0) and stepped (v). 



that reaches zero friction at a nonzero area. However, the 
friction is clearly not zero below this point and results at 
lower area show pronounced curvature. The curvature is 
even more pronounced when Oo is used instead of as 
would be expected from the fact that scales as load 
to the 1/3 power with an offset. We conclude that when 
a film is present the friction force is more strongly con- 
nected to normal load than area of contact. 



IV. SUMMARY AND CONCLUSIONS 

A single fluid monolayer was shown to have a sig- 
nificant impact on the mechanical properties of single- 
asperity contacts. Although the molecules diffused 
rapidly on an isolated substrate, they entered a glassy 
state when confined under the tip. As a result, the film 
was only expelled at the highest loads and for an in- 
registry commensurate tip where the barriers to motion 
were smallest. Previous simulations of hydrocarbons be- 
tween flat surfaces have also found films remain stable to 
very high pressures [63l [64] . The pressure gradient pro- 
duced by spherical tips should facilitate the expulsion of 
the film, but it remained stable over the times available to 
our simulations (~ 1/is). The loads used are comparable 
to the largest in AFM experiments and produced surface 
strains of more than 10%, which would lead to plastic in- 
dentation for most potentials other than the ideal springs 
considered here. 

Sliding at high load was more likely to damage the film, 
but the disruption rapidly healed through fluid diffusion 
once the tip had passed. Sliding only produced direct 
contact for the stepped tip and only at the highest load. 
The in-registry commensurate tip expelled the film under 
static loading, but lifted on top of the film during slid- 
ing. All the above results suggest that relatively weakly 



adsorbed molecules may remain within contacts under 
typical loading conditions. Their ability to separate op- 
posing surfaces and recover from damage are desirable 
features for boundary lubricants. 

The adsorbed film spreads the contact pressure over a 
substantially larger area (Figs. [5] and [7| than for bare tips 
or Hertz theory. Deviations in the variation of pressure 
with radius between different tip geometries were smaller 
than for bare tips because the film smeared features out 
over a molecular diameter. This smearing also explains 
the large increase in the outer radius of the pressure dis- 
tribution relative to Hertz. Values of are shifted by a 
constant of order the radius where the tip surface has 
risen above its lowest point by a molecular diameter, 
EE y/2aR. This is a significant fraction of the radius 
for typical AFM tips, but would be a small effect in the 
larger contacts of a Surface Force Apparatus. 

The pressure in the outer regions is small and Qq may 
overestimate the effective size of contacts. Contact radii 
derived from moments of the pressure distribution are 
closer to the values for bare tips and to Hertz theory. At 
large loads, is shifted from Hertz theory by a constant 
offset of a few molecular diameters. As N decreases, 
saturates. In this regime the local contact pressure is 
comparable to the ideal gas pressure. The tip is sup- 
ported by thermal fluctuations in the upper edge of the 
monolayer. The width of this layer is of order a and the 
pressure is spread over a radius of order Co- . 

The contact area was also determined from the num- 
ber of tip atoms in contact at any instant multiplied 
by the area per atom. The corresponding radius is 
just slightly below the Hertz prediction for most tips and 
loads. The source of the large discrepancy between 
and fls at low loads is that each tip atom is only in con- 
tact for a small fraction of the time. This is consistent 
with the tip being supported by thermal fluctuations in 
the upper edge of the monolayer. 

A surprisingly large temporal variation in the identity 
of the atoms that make contact was also observed at high 
loads. The number Nc of tip atoms that contacted the 
adsorbed film over a time interval t grew roughly as t^^'^ 
while the number contacting a bare substrate grew as In t. 
For bare surfaces the time dependence appears to reflect 
thermal vibrations in the substrate. For the adsorbed 
film, slow displacements of molecules to new glassy states 
and fluid fluctuations in low pressure regions at the edge 
of the contact also increase Nc- The time dependence 
of Nc makes it a more ambiguous measure of contact 
area than measures based on the mean pressure. Even 
increasing t to the period of the fastest phonon vibra- 
tions changed the contact area by 50%. Purely repul- 
sive interactions were used here, but an attractive tail in 
the potential also complicates the definition of contacting 
atoms. 

The adsorbed film also increases the normal tip dis- 
placement S, which is directly accessible in some exper- 
iments |12j . The in-registry commensurate tip exhibits 
a large jump in displacement when the film is expelled. 
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For other tips the displacement rises at the rate predicted 
by continuum theory at large loads. At smaller loads 5 
rises more rapidly. This is qualitatively consistent with 
continuum theory for a thin elastic coating where the ef- 
fective modulus crosses over from that of the coating to 
that of the more rigid substrate as the contact radius be- 
comes much larger than the film thickness [SHI ED El] ■ 
The deformation of the substrate is very close to Hertz 
theory, but is slightly reduced because the pressure is 
spread over a greater area. 

The normal displacement is closer to Hertz theory than 
the contact radius because it reflects an average response 
of the entire system, while the contact pressure is a local 
property. Fits at large loads give accurate values of the 
substrate modulus, indicating that analogous experimen- 
tal measurements provide accurate material parameters. 

The static and kinetic friction were studied for tips 
moved at constant velocity or pulled by a spring. The 
static friction is insensitive to the stiffness of the system, 
but the kinetic friction decreases as the stiffness increases 
(Fig. 111. Similar behavior is observed in the Prandtl- 



Tomlinson model |501 1541 I55j and previous simulations 
of tips moving over a fixed potential dHl ISSl ISZI- The 
rate of decrease is greatest for commensurate tips where 
atoms can move coherently between metastable states 
and the tip feels a nearly fixed potential. Amorphous 
tips show the least decrease because they can not jump 
to an equivalent minimum. The resulting rearrangements 
in the film lead to greater dissipation and a larger Fk- 

For all tips, the static friction is roughly proportional 
to N at low loads. Nonlinear behavior sets in at high 
loads when there is a transition, such as contact of a new 
terrace on a stepped tip or disruption of the film. The 
kinetic friction follows the same trends and is slightly 
smaller than Fs when the tip is pulled through a weak 
spring. For a stiff system, the friction is extremely small 
for a range of loads and then rises linearly. The extent of 
the low friction region is largest for commensurate sys- 
tems because of the coherent motion discussed above. 
The value of Fk in this region goes to zero with decreasing 
velocity. A similar low friction region has been observed 
in AFM experiments at low loads and is consistent with 
the Prandtl-Tomlinson model [501 EH [SS] . The apparent 
discrepancy between previous results for friction on bare 
amorphous tips [121 EZl HE] could be explained by the 
fact that one reported Fg while the latter reported F^. 

The relation between friction and contact area was also 
examined for each of the definitions described above. The 
best correlation was obtained for the static friction and 
(Fig. 13), but the dependence is less linear than for 



load, and appears to go to zero at finite area. Note that 
scales roughly as N"^/^ (Fig. [?]). If friction is linear 
in load, plots of friction against area should have the 
upwards curvature observed in Fig. |13[ 

The linear relation between friction and load is con- 
sistent with previous simulations of adsorbed layers be- 
tween flat surfaces. For a wide range of atomic structure 
and surface coverage the shear stress Tshear is a linear 
function of the local pressure: Tgj^^ar = ''o + otp. If this 
linear relation persists between rough surfaces, then the 
total friction should scale as F = T^Areai + otN for any 
distribution of pressure. For the case of nonadhesive in- 
teractions considered here tq is nearly zero and friction 
should rise linearly with load as observed. Moreover a di- 
rect comparison of the values of a for simulations of flat 
surfaces with the ratios Fg/N observed here shows the 
same trends with structure and similar numerical values. 
For example, commensurate surfaces yield a ~ 0.3, in- 
commensurate surfaces have a « 0.05, and amorphous 
surfaces have a ~ 0.1 [251|271I2S] Simulations also show 
that a is insensitive to the exact coverage |29| . explain- 
ing why damage to the film does not produce dramatic 
changes in F/N. 

Studies of thin films between surfaces with a small ran- 
dom roughness found a similar linear rise in Tshear |47j . 
The authors argued that the real area of contact was ir- 
relevant in these simulations because there was no direct 
contact of opposing surfaces, but any definition based on 
contact pressure would give the real area of load propor- 
tional to the apparent area in their simulations. It would 
be interesting to extend simulations to much rougher sur- 
faces where the contact pressure is only significant over a 
small fraction of the apparent area and where this Areai 
rises with load. 

Another important extension will be to consider a 
monolayer which partially wets the tip. In this case a 
capillary meniscus will form around the tip. In addition 
to producing an adhesive force on the tip, the capillary 
may introduce new dissipation mechanisms. It may also 
serve as a reservoir that feeds molecules to the contact 
and prevents film rupture. Another parameter that could 
be varied is the thickness of the film on the substrate. 
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